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ABSTRACT 

We derive a physical model for the observed relations between star formation rate (SFR) and molecular line 
(CO and HCN) emission in galaxies, and show how these observed relations are reflective of the underlying 
star formation law. We do this by combining 3D non-LTE radiative transfer calculations with hydrodynamic 
simulations of isolated disk galaxies and galaxy mergers. We demonstrate that the observed SFR-molecular 
line relations are driven by the relationship between molecular line emission and gas density, and anchored by 
the index of the underlying Schmidt law controlling the SFR in the galaxy. Lines with low critical densities 
(e.g. CO J=l-0) are typically thermalized and trace the gas density faithfully. In these cases, the SFR will be 
related to line luminosity with an index similar to the Schmidt law index. Lines with high critical densities 
greater than the mean density of most of the emitting clouds in a galaxy (e.g. CO J=3-2, HCN J=l-0) will have 
only a small amount of thermalized gas, and consequently a superlinear relationship between molecular line 
luminosity (Lmoi) and mean gas density («). This results in a SFR-line luminosity index less than the Schmidt 
index for high critical density tracers. One observational consequence of this is a significant redistribution of 
light from the small pockets of dense, thermalized gas to diffuse gas along the line of sight, and prodigious 
emission from subthermally excited gas. At the highest star formation rates, the SFR-Lmoi slope tends to 
the Schmidt index, regardless of the molecular transition. The fundamental relation is the Kennicutt-Schmidt 
law, rather than the relation between SFR and molecular line luminosity. Our model for SFR-molecular line 
relations quantitatively reproduces the slopes of the observed SFR-CO (J=l-0), CO (J=3-2) and HCN (J=l-0) 
relations when a Schmidt law with index of ^ L5 describes the SFR. We use these results to make imminently 
testable predictions for the SFR-molecular line relations of unobserved transitions. 

Subject headings: galaxies: ISM - galaxies: starburst - ISM: molecules - stars: formation - radio lines: 
galaxies - radio lines: ISM 



1. INTRODUCTION 

The rate at which stars form in galaxies has historically 
been parameterized in terms of "laws" relatin g the star forma- 
tion rate (SFR) to the density of available gas. lSchmidJ (Il959l) 
originally proposed a power-law form for the SFR such that 
SFR oc p'^ (hereafter, referred to as a Schmidt Law). 

Observed SFR relations typically come in two flavors. The 
first, relating surface SFR density to surface gas density takes 
the form: 



Esfr oc !]„ 



(1) 



Observations of local galaxies have constrained t he surface 
density SF R index, N, to roughly A^= L4±0. 1 5 (e.g. lKennicuttI 
[l998a b: Kennicutt et al.ll2007. and references therein). 

The second varietal of SFR indicators relates the SFR to 
the mass of molecular gas above a given volumetric density. 
For example, local galaxy surveys have shown a relation ex- 
ists between the SFR and molecular gas such that the SFR 
(as traced by the infrared luminosity - hereafter Lir) is pro- 
portional to the ^^CO (J=l-0) luminosity to the 1 .4- L6 power 
dSanders et al.lll991l:ISanders & Mirabell 19961 and references 
therein). Because the J=l-0 transition of '^CO (hereafter CO) 
can be excited at relatively low densities (n ^10^-10"' cm"^), 
and lies a modest ^5 K above ground, it serves as a valuable 
tracer of total molecular H2 gas content down to relatively 
low densities. These observed relations have been broadly in- 
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terpreted as an increasing star formation efficiency (SFE; the 
SFR divided by Mm) a s a function of molec ular gas mass 
(iGao & Solomonl2004allbtlSanders et alWmi S. 

Recent observations of HCN have suggested that a perhaps 
more fundamental volumetric SFR relation for galaxies exists 
in terms of the dense molecular gas. Because HCN (J=l-0) 
has a relatively high critical density («crit ^10^ cm~^) com- 
pared to that of CO (J=l-0), HCN is typically only thermal- 
ized in the dense cores of molecular clouds. Thus, in the limit 
that the bulk of HCN luminosity originates from thermalized 
gas, HCN serves as a good tracer of the dense molecular gas 
that is actively involved in the star formation process. This 
is in contrast to CO (J=l-0) which tends to emit from both 
dense cores as well as diffuse molecular filaments and cloud 
atmospheres. 

Pioneering millimeter-wave observations by 

iGao & SolomonI (l2004albl) uncovered a tight linear cor- 
relation between the infrared luminosity and HCN luminosity 
in a sample of local galaxies, ranging from quiescent spirals 
to luminous and ultraluminous infrared galaxies (LIRGs 
an d ULIRG s , resp ectively). This relationship was found 
by IWu et all (l2005h to hold in individual star forming cores 
in the Milky Way as well. Observations of a roughly 
linear correlation between Lir and CO (J=3-2) emission 
(with critical density n^it ^IQ'^ cm"-') in a similar sample 
of galaxies further corroborated this result, and provided 
evidence against HCN-re lated chemis try driving the observed 
Gao& SolomonI (l2004allbb relations dNaravanan et all 120051: 



Yao et alJl2003r 



Despite a plethora of observations, a consensus physical 
interpretation of the seemingly disparate Lir-HCN (J=l-0), 
LiR-CO (J=l-0) and Lir-CO (J=3-2) relations has yet to be 
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borne out of the literature. i Gao & SolomonI (l2004al[bl) sug- 
gest that the Hnear relationship between Lir-HCN (J=1-0) re- 
flects a scenario in which the SFR in galaxies is directly pro- 
portional to the amount of dense molecular gas, and simi- 
larly, a constant star formation efficiency in terms of dense 
molecular gas. In this picture, LIRGs and ULIRGs simply 
have a higher fraction of their molecular gas in a dense gas 
phase which is reflected in the linear Lir-HCN (J=1-0) slope^ 
and non-linear LiR-CO (J=l-0) slope. [Narayanan et al.l(l2005h 
suggest a similar interpretation for their observed relation of 
L trocCO(J=3-2)" '^^. 

IWu et al.l (l2005h extended these interpretations to include 
the observed linear Lir-HCN (J= 1-0) relation within the 
Galaxy. These observations suggested that dense molecular 
cores may represent a fundamental unit of star formation. In 
this case, the only difference between a ULIRG like Arp 220 
and a Galactic star forming region is the number of individ- 
ual star forming units present. This leads to a natural linear 
relationship between the SFR and amount of dense molecular 
gas. 



From the theoretical side, iKrumholz & ThompsonI (l2007h 
made headway towards understanding the observed relations 
by providing a model which quantitatively reproduces the ob- 
served LiR-CO (J=l-0) and Lir-HCN (J=1-0) relations for star 
forming clouds. By combining ID non-local thermodynamic 
equilibrium (LTE) radiative transfer calculations with models 
of giant molecular clouds (GMCs), these authors found that 
the relationship between the SFR and molecular line luminos- 
ity in star forming clouds depends on how the critical density 
of the molecule compares to the mean density of the observed 
source. Lines with critical densities which are typically below 
the mean density in a galaxy (e.g. CO J=l-0) probe the total 
molecular gas commonly from galaxy to galaxy. In this sce- 
nario, the SFR-molecular line exponent can be represented by 
the index that relates SFR to total volumetric gas density (i.e. 
the Schmidt index). Conversely, observations of molecules 
which typically trace densities well above the mean density 
of the galaxy (e.g. HCN J=l-0) trace similar conditions from 
galaxy to galaxy - i.e. the peaks in the density spectrum. In 
these cases, the molecular line luminosity rises faster than lin- 
early with increasing gas density and the corresponding rela- 
tion between SFR and line luminosity is near linear. A key 
direction forward beyond these models is understanding the 
effects of emitting GMCs on a galaxy-wide scale, with the 
potential effects of molecular line radiative transfer 

In this paper, we build on the body of observational and the- 
oretical work by providing a physical model for the origin of 
the observed molecular SFR indicators on galaxy-wide scales, 
and relating them to observed relations. We additionally make 
model-distinguishing testable predictions for how the SFR in 
galaxies relates to as yet unobserved transitions in CO and 
HCN. We do this by combini ng 3D non -LTE molecular line 
radiative transfer codes ( Naravanan et al...2006a^K .2007 ab.) 
with smoothed particle hydrodynamic (SPH) simulations of 
both isolated star forming galaxies and galaxy mergers. Our 
methodology includes the effects of both collisional and ra- 
diative molecular excitation and de-excitation, a m ulti-phase 
ISM, and star formation (IS pringel. Di Matteo & Hernquist 
l2005h . We additionally include a methodology for black hole 
growth and the winds associated with both black hole growth 

^ Throughout this paper we use the term slope interchangeably with 'in- 
dex' (as in the exponent in a given SFR relation - e.g. Equation [T). We do 
this as we are considering all relations in log-log space. 



and star formation, though note that these processes play a 
minimal role on the SFR-L,noi relation in our simulations. 

The paper is outlined as follows. In § |2] we describe our 
hydrodynamic and radiative transfer simulations. In §|3]we 
quantitatively describe the origin for the observed SFR-CO 
and SFR-HCN slopes and follow by making testable predic- 
tions for submillimeter-wave telescopes in § |4] In § |5] we 
compare our results to the body of observational work in this 
field, and in § |6] relate these simulations to other models and 
interpretations for molecular SFR indicators. In §|7] we con- 
clude with a summary. Throughout the work we assume a 
ACDM cosmology with /i=0.7, r2A=0.7, fiM=0.3. 

2. NUMERICAL METHODS 
2.1. Hydrodynamics 

For this work we have modeled both iso lated disk galax- 
ies and major mergers. We do this a s the iGao & SolomonI 
(2004a'.D) and 'Narav anan et al.l (1200 5') samples include both 
disk galaxies as well as ongoing mergers. The merger snap- 
shots are spaced equally temporally with a sampling of 5 Myr. 
We note, however, that the results in this paper are not de- 
pendent on our usage of any particular combination of disk 
galaxies or mergers. 

The hydrodynamic simulations were conducted with a 
modified version o f the pubhcly available A^-body/SPH 
code GADGET-2 (ISpringell l2005h . The prescriptions 
used to generate the galaxies, as well as the algo- 
rithms involved in simulating the physics of the multi- 
phase ISM, star formatio n, and black hole growth a re de- 
scri bed fuUy in in IS pringel. Di Matteo & Hernquist! (l2005l) 
and ISpringel & Hernquist (.2002. .2003.) . We refer the reader 
to these works for further details though summarize the as- 
pects most relevant to this study here. 

GADGET-2 accounts for radia tive cooling of the gas 
(lDaveetal.lll99l iKatz et al.1 [19961) . and a multi-phase ISM 
which is considered to consist of cold clouds embedded i n 
a hot, pressure confining ISM (e.g. lMcKee & Ostrikedll9'77[) . 
This is realized numerically through 'hybrid' SPH particles in 
which cold clouds are allowed to grow through radiative cool- 
ing of the hot ISM, and conversely star formation can evapo- 
rate the cold medium into diffuse, hot gas. Pressure feedback 
from supernovae heating is t reated via an effective equation of 
state (EOS) (see Figure 4 of'Spring el. Di Matteo & Hernquist! 
!2005h . Here we utilize an EOS softening parameter of 
q'EOS=0.25. 

Star formati on is constrained to fit the observe d Kennicutt- 
Schmidt laws (!Kennicutt!!l998a![bt!SchrnId3!l959l) . In particu- 
lar, star formation is assumed to be related to the gas density 
and local star formation time scale: 
dpi. 



dt 



K1-/3) 



Pcold 



(2) 



where (3 quantifies the mass fraction of stars which are short- 
lived and effectively immediatley supernova, and f* is the 
timescale for star formation. 

A Salpeter IMF with slope -L35 and mass range 0.1-40 
M0 returns a value of (3 ^0. 1 . The parameter f* is assumed 
to be proportional to the local dynamical time scale. 



t*ip) = t{) 



-1/2 



(3) 



Detailed sfiidies by !Cox et al.! (!2006cl) and 
!Springel & Hernquist! (!2003h have found that a propor- 
tionality constant of fj = 2.1 Gyr reproduces both the 
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normalization and slope of the iKennicutl (Il998al[bl) observed 
SFR surface density relations over mergers with a wide 
variety of isolated galaxy and galaxy merger models. The 
dispersion in the modeled relations lay well within the 
observed dispersion . Simi lar resu lts have been fo und by 
iMihos & HernquistI (ll994allR IT996) and Springell (12000) . 
Hence, we use this formulation for the star formation rate. 

Black holes are optionally included in our simulations 
which accrete via a Bondi-Lyttleton-Hoyle parameterization 
with a fixed maximum rate corresponding to the Eddington 
limit. The black hole radiates such that its bolometric lumi- 
nosity is set by the accretion rate with L = eMc^ with accretion 
efficiency e=0.1. We further assume that 5% of this energy 
coupl es isotropically to the surrounding ISM as feedback en- 
ergv (iDi Matteo et al.ll2005l:ISpringel. Pi Matteo & HernquistI 
12005b . a number chosen to match the normalization of the 
locally observed MsH-Mbuhe relation. We have run identical 
simulations both with and without black holes for this work, 
and found that the inclusion of black holes makes little differ- 
ence on the final result. We note, however, that we focus our 
molecular line comparisons to the SFR, and not Lir (which 
could, in principle, have a contribution from buried AGNs at 
the highest luminosities) as a complete treatment o f model- 
ing the infrared SEP (e.g. lChakrabarti et al.ll2007allbl : lJonssonl 
l2006t iLi et al.ll2007bh is outside the scope of this study. 

We consider both isolated disk galaxies with varying gas 
fractions and masses, as well as two equal mass binary merg- 
ers. We utilized 15 isolated galaxies exploring a parame- 
ter space with gas fractions=[0.2, 0.4, 0.8] and virial veloc- 
ities=[115, 160, 225, 320, 500 kms"']. All disk galaxies 
(including the pr ogenitors for the merger simulations) were 
initialized with a iHernquislI d 19901) dark matter profile, con- 
centration index c=9, and spin parameter A=0.033. For all 
galaxies we utilize 120,000 dark matter particles, and 80,000 
total disk particles. The softening lengths were 100 pc for 
baryons and 200 pc for dark matter 

For the mergers, the progenitor galaxies were identical to 
the disk galaxies described here, with V2oo=160 km s"', and 
40% gas fraction. For reference, the mer gers are simula- 
tions 'n o- winds' and 'BH' as described in i Naravanan et all 
(l2007bl) . The first merger snapshots considered are when 
the molecular disks of the progenitors are overlapped as the 
galaxies approach final coalescence. We utilize multiple snap- 
shots throughout the evolution of the merger simulations be- 
yond this point to represent observed merging pairs caught in 
various stages of evolution (e.g. close pairs such as the An- 
tennae to more evolved mergers such as Arp 220 and NGC 
6240). All merger snapshots are during the period of high- 
est star formation activity, when the galaxy can be consid- 
ered a LIRG/ULIRG (model IR lumino sities for these merger 
simulations can be found in Chakrabarti et al.l (l2007al) ). Sam- 
ple CO images of the merger simulations employed here, as 
well as some of the m odel disk galaxies can be found in 
INaravanan et all (l2007bl) . 

It is an important point to note that no particular choice or 
combination of models affects the presented results. While 
some of the disk galaxies used may nominally have larger cir- 
cular velocities or gas fractions than galaxies observed in the 
local Universe, because the results are general, not including 
them does not change our results or interpretation. We include 
a wide range of model galaxies to increase the number statis- 
tics and dynamic range of our simulations. The generality of 
our results will be shown more explicitly in §|3]and §|4| 



2.2. Non-LTE Radiative Transfer 

The propagation of a line through a medium with lower 
density than the line's critical density requires a full non-LTE 
treatment. Specifically, in this regime, both collisions and ra- 
diative processes contribute to the excitation and de-excitation 
of molecules. This formalism for molecular line transfer has 
long been applied to interpretations of galaxy observations 
via la rge velocity gradient codes (e.g. 'Goldrei ch & KwanI 
1 19741) . though has only recently b een incorp orated in full 3 
dimensions in galaxy-\vide models (|Greye & S ommer-Larsei 
2006t INaravanan et al.' '2006al l2007allbl: IWada & Tomisak; 



2005 



; [Yamada et al. 2007 ). We employ the radiative transfer 



methodology of Naravana n et al.l (l2007b). and refer the reader 
to that work for full details. Here, we briefly summarize. 

The radiative transfer is performed in two phases. First, the 
molecular level populations (and consequently source func- 
tions) are ex phcitly calculated using the 3P non-LTE code, 
Turtlebeach (INaravanan et al.l l2006a' '2007b'). In this, a solu- 
tion grid to the level populations is guessed at, and model pho- 
tons are emitted isotropically in a Monte Carlo manner The 
level populations are updated by assuming statistical equilib- 
rium, and balancing radiative and c ollisional exc itations, de- 
excitation and stimulated emission (lBernesl[l979l) . A new se- 
ries of model photons is emitted and the process is repeated 
until the level populations are converged. 

Once the level populations and source functions are known, 
the model intensity can be found by integrating the equation 
of radiative transfer along various lines of sight. Formally: 



ro 



SAr)[l- 



-r„(r)l 



-T„ (tot) 



(4) 



where I^ is the frequency-dependent intensity, S^, is the source 
function, r is the physical depth along the line of sight, and r 
is the optical depth. 

The SPH outputs are smoothed onto a grid with ^^160 pc 
spatial resolution for the radiative transfer. In order to more 
accurately describe the strongly density-dependent collisional 
excitation and de-excitation rates, we model the gas in grid 
cells to be bound in a mass spectrum of GMCs constrained 
by observations of Milky Way clouds (Blitz et al. 2006). The 
GMCs are modeled as spheres with power-law density gradi- 
ents, and r adii given by the Galactic GMC mass-radius re - 
lation (e.g. iRosolowskvl l200l 120071: ISolomon et all [19871) . 
These GMCs are placed randomly in the grid cells (which we 
refer to as "cells" hereafter). Of the SPH simulations utilized 
here in particular, the maximum volume-averaged H2 density 
seen in a cell was of order ^1000 cm"^, though higher max- 
imum densities are of course realized in the subgrid GMCs, 
dependent on the mass of the GMC (typically central core 
densities reached ^10''- 10^ cm"^ in the nuclear regions of star 
forming galaxies in our simulations: iNaravanan et alj|2007bl) . 
Because the sub-grid GMCs are posited on the grid in post- 
processing, the SFR (Equation|2l) calculated during the hydro- 
dynamic simulations utilizes the volume-averaged densities 
achieved at the hydrodynamic resolution (which is of order ^ 
100 pc). More details concerning the subgrid formalism for 
including GMCs on a galaxy- wide scale, and resolution tests 
may be found in ( INaravanan et alj|2007bl) . 

Observational evidence suggests a range of power-law in- 
dices for GMCs, ran ging fr om n=\-2 (Andre etal. 199^ 
iFuller & Mverslll992i: IWalker. Adams & Ladal ll990). Tests 
utilizing a number of cloud density power-law indices within 
this range (in which the total mass of the cell was con- 
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served, and thus the central density in the clouds allowed to 
vary) show that the results in this work are not sensitive to 
this param eter choice. Similarly, in the study of iBlitz et alj 
(12006) and Rosolowskvl (12007). GMC mass spectrum indices 
were found to range from 7 «-1.4 to -2.8. Again, tests of 
mass spectrum indices showed that the results of this pa- 
per are not sensitive to parameter choices within the range 
of observational constraints. We nominally employ indices 
of n=1.5 and 7=-!. 8 for the cloud dens ity and GMC mass 
spectrum power-laws ( lRosolowskvll2007h . This methodology 
allows us to faithfully capture the emission processes from 
both dense molecular co res and more diffuse GMC envelopes 
(iNarayanan et al.ll2007bh . 

We have benchmarked our radiative transfer codes 
against published non- LTE radiative transfer tests 
(Ivan Zadelhoff et aP l2002h, and pr esent the results for 
these tests in lNaravanan et al] (l2006bh . Our methodology for 
applying 3D non-LTE molecular line radiative transfer has 
shown success in reproducing characteristic observed CO line 
widths, morphologies, excitation conditions and intensities 
in isolat ed disk galaxies, loc al ULIRGs and quasars from 
z ~ 2-6 dNaravanan et al.ll2006 a, 2007a b). 

The molecular gas mass fraction is a ssumed to be half, as 
motiv ated by local volume surveys (e.g. lKeres. Yun & Yound 
l2003h . and the molecular abundances set uniformly at Galac- 
tic values (Lee, Bettens & Herbst 1996). In this work, typ- 
ically '^IxlO^ model photons were emitted in the non- 
LTE calculations per iteration, and we considered transitions 
across 1 1 molecular levels at a time. The mass spectrum of 
GMCs had a lower mass cutoff of 1 x lO'^Mo and an upper 
mass cutoff of 1 x lO^M© .. The collisional rate coefficients 
were taken from the Leiden Atomic and Molecular Database 
(ISchoier et al.ll2005h . 

3. ORIGIN OF OBSERVED SFR-CO AND SFR-HCN SLOPES 

3.L General Argument 

In Figure [T] we plot the SFR-CO relations and SFR- 
HCN (J=l-0) relation as derived from our models of iso- 
lated disk galaxies and mergers with best fitting slopes over- 
laid. Here, we show the emission from a sample of model 
galaxies as a whole (e.g. simulated unresolved observa- 
tions of galaxies). We plot the line intensities, though note 
that they are proportiona l to m olecular line luminosity (L', 
ISolomon & Vanden Bout! (l2005l) ). the standard quantity re- 
ported in the observational literature*. Included in the plot 
are a random sampling of 35 of the 15 disk galaxies and ^80 
merger snapshots in our parent sample. The points are all 
different simulated galaxies, and not simply multiple viewing 
angles of the same galaxy. We first note the general agreement 
of the best fit slope s and disp e rsion with observed rel ations 
by iGao & Solomon' (l2004albi) . iNaravanan etall (12005 ^. and 
[Sanders et al. ( 1991). The dispersion in Lmoi at a fixed SFR 
arises owing to variance in the mean density of the emitting 
gas in galaxies at a given SFR. We will discuss the physi- 
cal reasoning behind the slopes shortly. Second, we note that 
we have utilized a randomly drawn sample of our simulated 
galaxies to demonstrate the generality of our results. From 
this, we see that any given combination of model disk galax- 
ies and/or merger snapshots give roughly the same fits. We 
have quoted as "error" values the \a dispersion in 100 ran- 

' L' Ri 7r/(4 In 2)d^_^^I^^al'^(^+zy^ where 9„,5 is the main beam angular 
size, /,noi is th e molecular line intensit y, dt is the luminosity distance, and z 
is the redshift IGao & Solomonl2004al) . 



dom draws of 35 model snapshots. The generality of these 
results occurs because the origin of the SFR-molecular line 
relations arise from the nature of the molecular emission from 
the individual galaxies themselves. 
To see this, consider a galaxy which is forming stars at rate: 



SFR ex / 



(5) 



where p is the mean molecular gas mass density. The rela- 
tionship between the SFR and the luminosity of a molecular 
line from a cell of clouds 



SFR oc l: 



'molecule 



(6) 



is dependent on the relationship between the molecular line 
luminosity and the mean gas density in a cell of clouds; 



^molecule ^ P 



(7) 



where (3=N/a. Therefore, for a given Schmidt law index, A^, 
the root issue in determining how the SFR relates to molecular 
line luminosity is understanding how the molecular line traces 
molecular gas of different densities (/3). 

The global relationship between molecular line luminosity 
and mean gas density traced (/3) is driven by how the critical 
density of the molecular line compares with the density of the 
bulk of the clouds across the galaxy. In short, lines which 
have critical density below the mean density of most of the 
emitting gas cells will be thermalized and rise linearly with 
increasing cloud density. These lines will consequently have 
SFR-line luminosity relations (a) similar to the Schmidt index 
controlling the SFR. 

Lines which have critical density well above the mean den- 
sity of most of the emitting gas cells will be thermalized in 
only a small fraction of the gas. To see this, in Figure |2] 
we plot the distribution of mean cell densities for a fiducial 
disk galaxy (with 40% gas fraction and 160 km s"' circular 
velocity; this will be a fiducial disk galaxy which we center 
the remainder of the discussion around as a reference point, 
though the results are general)^. We additionally plot the crit- 
ical densities of a sample of CO lines overlaid as an indicator 
of the relative gas fractions above given critical densities. In 
the case of high critical density tracers such as e.g. CO (J=3- 
2), a small fraction of the gas is thermalized. As the mean 
density increases, the fraction of thermalized gas (and con- 
sequently photon production) increases superlinearly, driving 
an e.g. SFR-CO (J=3-2) relation with index less than that 
of the Schmidt index. This is similar to the physical mecha- 
nism driving the SFR-L,noi relation in star forming GMCs by 
iKrumholz & Tho mpson (2007). 

An important question to understand is, how does this 
global model for the observed SFR-L r rmi rel ations (also con- 
sidered by IKrumholz & Thompson! |2007|) manifest itself 
when considering the actual transfer of molecular line photons 
in physical models of isolated galaxies and galaxy mergers? 
In particular, high critical density photons from the dense, 
thermalized nucleus may be reabsorbed by diffuse gas along 
the line of sight. How, then, considering this redistribution, is 
the observed SFR-Lmoi relation realized? In short, in galax- 
ies where n ^ «ciit, large amounts of diffuse gas along the 
line of sight are subthermally (radiatively) excited owing to 
line trapping. Emission from subthermally excited gas can be 

^ In order to account for the destruction of molecules in photodissociation 
regions (e.g. HoUenbach & Tielens 1999), we do not consider emission from 
regions with column density ^l-5x 10-' cm"-. This typically corresponds to 
cells with mean cloud density ~50 cm"' in our simulations. 



Molecular SFR Indicators in Galaxies 





CO 1-0 








' 




.. / 


100 




• 


•V ■ 






• / 


/ ; 




/ 


/ • 


• 


10 


/ 


• 





C03-2 



HCNl-0 



100 



10 



10-" 



10'" 
Ico 



10' 



10" 




100 



10 



SFR ocI„ 



10-" 



10-" 



10-'^ 



10-' 



Fig. 1 . — Model results showing the relationship between SFR and CO (J=l-0), CO (J=3-2) and HCN (J=l-0) emission in simulated galaxies. The molecular 
line emission is derived for a randomly drawn setof ~35 galaxy snapshots of our model sample of ~ 100, and is a velocity-integrated line intensity. The randomly 
drawn set is different for each plotted transition in order to illustrate the generality of our results. The solid lines are the least-squares fit to the plotted model 
results, and the fit is quoted in each panel. We randomly drew the sample of 35 galaxies 100 times, and quote the Icr dispers ion in the derived slopes as the 
'error' values. The modeled relationships between SFR and molecular line luminosity are consistent with the observations of Gao & Solomon' l2004a b) and 
[Narayanan et al. (2005). This figure takes the emi ssion from t he simulated galaxies as a whole, thus simulating unresolved detections of the galaxies. Thus it 
serves as a direct comparison to the observations bv lGao & Solomon (i2004ai ibi) and iNaravanan et alj)2005) . While this particular plot shows simulated unresolved 
observations of galaxies, we refrain from using traditional molecular line "luminosity" units as future discussion will center largely around the physical intensity 
associated with individual emitting gas cells in the simulations. The units of intensity are in erg s"' cm"^ Hz"' . 
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Fig. 2. — Distribution of mean cloud densities throughout our fiducial disk 
galaxy with circular velocity i'=160 km s"' and gas fraction /=0.4. The hues 
overlaid denote the location of critical densities of common CO transitions, 
and the color of each line indicates the fraction of gas mass with density 
above each critical density (with color scale on right). While most of the 
gas mass is above the critical density of the lowest lying transition, CO (J=l- 
0), higher lying transitions (e.g. CO J=3-2) may be subthermally populated 
through much of the galaxy. 

thought of as a scattering process in that the high critical den- 
sity photon originates in dense, thermalized molecular gas, 
and is redistributed to more diffuse gas. Subthermal emis- 
sion from radiatively excited cells tends to have a superlinear 
Lmoi-n relation as the relatively diffuse gas owes its luminos- 
ity to molecules who have been pumped by radiation from 
denser regions. As such, the radiation is indicative of denser, 
thermalized regions where the photons were originally cre- 
ated, but the diffuse gas along the line of sight reduces the 
sightline-averaged mean density. In these cases, the SFR-line 
luminosity relation (a) has an index less than that of the un- 
derlying Schmidt law for the galaxy. There are details, of 
course, specific to the relationship between CO (J=l-0), CO 
(J=3-2) and HCN (J=l-0) luminosities and gas density traced 
which require a more in depth analysis, but these underly- 
ing themes are robust. This is inherently an alternative but 
equivalent way of viewing the previous point. The net num- 
ber of high critical density photons depends on the quantity of 
dense, thermalized gas. The bulk of the gas along the line of 



sight (which is responsible for the line trapping) is diffuse, re- 
sulting in a superlinear Lmoi-n relation along the line of sight. 
We devote the remainder of this section to examining the 
driving mechanisms behind the relationship between line lu- 
minosity and mean gas density (/?) in greater detail, and uti- 
lize these to formulate a general model for observed SFR- 
molecular line luminosity relations in galaxies. We will utilize 
the aforementioned fiducial disk galaxy as a reference point. 
We will couch the discussion (at least for the high critical den- 
sity tracers - e.g. HCN J=l-0 and CO J=3-2) in terms of both 
the global relationship between h and «crit, as well as the redis- 
tribution of photons from thermalized gas to diffuse, subther- 
mally excited cells. The reason we do this is the following. 
The SFR-Lmoi relation in galaxies is globally set by the num- 
ber of high critical density photons produced as a function of 
the mean density of the emitting clouds along the line of sight 
(e.g. the relationship between h and ncrit; Figure |2]i. While 
these high critical density photons are originally produced in 
the small number of thermalized cores in galaxies, this does 
not necessarily reflect the final emitting surface in the galaxy. 
The photons are redistributed from dense cores to being ab- 
sorbed and re-emitted by diffuse gas along the line of sight. 
The details of how the SFR-L^oi relation is driven in a real- 
istic model for a galaxy are subtle, and require more analysis 
than the global arguments presented thusfar Moreover, the 
large quantities of subthermally excited gas is a direct obser- 
vational consequence of this physical model for the observed 
SFR-line luminosity relations. For these reasons, we continue 
forward framing the understanding of the observed relations 
between SFR and line luminosity in terms of both global prop- 
erties as well as the subthermally excited gas along the line 
of sight. That said, the details outlined in the remainder of 
this section are just that. The remainder of the predictions 
presented in this paper are not dependent on the finer points 
regarding the excitation patterns in galaxies. 

3.2. Detailed Understanding of an Individual Galaxy 

We first note that a significant portion of the intrinsic lumi- 
nosity in the galaxy escapes. While the optical depths can be- 
come large locally, velocity gradients across the galaxy serve 
to shift native absorption profiles out of resonance with the 
emission line on larger scales. Specifically, we see ^60% of 
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the emission across most CO and HCN transitions escape the 
galaxy. The gas cells contributing to emission that escapes are 
of a fixed size, and distributed across the galaxy. These cells 
exhibit a wide range of mean cloud densities. We therefore 
focus on the properties of individual gas cells in our simula- 
tions (of all masses and densities) in order to build an under- 
standing of how their summed properties drive the global line 
luminosities-gas density relation, and the consequent SFR- 
molecular line luminosity relation. 

The observ ed SFR-molecular line relations 
JGao & SolomonI l2004a b: Naravanan et alj [2005) are 
representative of three regimes: one in which the mean 
density of the bulk of the emitting clouds is higher than 
the line's critical density; one in which the mean density is 
lower than the line's critical density; and one in which the 
mean density is lower than the line's critical density and 
the line considered is a ground-state transition. We explore 
these three cases here with respect to our fiducial disk galaxy 
(though note again that the results are general). 

In Figure |3] we show empirically how the line of sight 
velocity-integrated intensity from CO (J=l-0), CO (J=3-2) 
and HCN (J=l-0) relate to the sightline-averaged mean den- 
sity along 50^ sightlines for the fiducial disk galaxy^. The 
solid line in each case shows linearity, with arbitrary normal- 
ization. In Figure m we plot the line luminosity versus mean 
density on a cell by cell basis for the CO (J=l-0), (J=3-2) and 
HCN (J=l-0) transitions from our fiducial disk galaxy. We re- 
mind the reader that a 'cell' is a grid cell in our fiducial galaxy 
which contains potentially num erous GMCs whose effec ts are 
simulated in a subgrid manner (Naraya nan et al.l l2007b). We 
will refer to these two figures throughout the forthcoming dis- 
cussion. 



3.2.1. iici-it <n 

We first consider the regime in which the mean density 
of the bulk of the emitting clouds is typically higher than 
the line's critical density. We see empirically from the left 
panel of Figure |3]that the CO (J=l-0) luminosity traces sight- 
lines of increasing mean density linearly in our fiducial galaxy 
if] = 1). Utilizing equations |5]|7] this translates into a globally 
observed relation in a sample of galaxies in which the SFR re- 
lates to Lmoi with an index similar to the underlying Schmidt 
index (here, A^=1.5). What causes this linear relation between 
Lmoi and n? 

At the critical density, the line is approaching thermaUza- 
tion and, in the limit that constant molecular abundances ap- 
ply (with respect to H2), the clouds are typically optically 
thick. In this scenario, the emission arises primarily from the 
outer (lower column density) regions of clouds, and serves 
as a measure of the number of emitting clouds in a given 
cell. While groups of clouds are locally optically thick, on 
a galaxy-wide scale they are optically thin. Thus, the veloc- 
ity dispersion of the bulk of clouds along a given line of sight 
has the effect of spreading the molecular line emission out in 
frequency space, permitting physically overlapping clouds to 
be counted. Because the mean cloud density increases with 
number of clouds in a given cell, in the optically thick limit, 
a roughly linear relation between molecular line emission and 
mean cloud density is natural. 

** Note that this is now plotting the L„,oi-n relation for individual sight- 
lines for a single disk galaxy (with implied spatial resolution equal to that 
of the simulation of ~160 pc). This is in contrast to Figure [T] which plotted 
simulated unresolved observations of galaxies as a whole. 



In Figure |4] (left panel), we illustrate this effect by showing 
the CO (J=l-0) intensity as a function of mean cell density in 
our fiducial galaxy. The points are individual grid cells con- 
taining (potentially) numerous GMCs. Overlaid in color is the 
level of thermalization. The bulk of the CO (J=l-0) emission 
comes from cells that are either in LTE or nearly thermalized. 
We therefore arrive at the conclusion that transitions that have 
critical density well below the mean critical density of most 
emitting clouds will be roughly thermalized, and the SFR-line 
index will be similar to the underlying Schmidt-law index. 
When an underlying Schmidt index of A^=1.5 is in place, the 
SFR-Ln,oi relation matches obs ervations well. This is similar 
to the explanation posited by iKrumholz & Thompson! (l2007h 
for their models of GMCs. 

). 

3.2.2. nci-it>fi 

If the critical density of the transition is toward the high 
density tail of the gas density distribution function, as one 
considers galaxies of increasing mean density, the luminos- 
ity from that molecular transition will increase superlinearly. 
This owes to an increase in thermalized gas both from an in- 
creased quantity of gas in the galaxy, as well as an increased 
fraction of gas above the critical density. This effect is seen in 
simulated global unresolved observations of galaxies (which 
show effectively a single, unresolved sightline in a galaxy) as 
well as simulated higher resolution observations. We illus- 
trate this intuitive argument in the middle and right panels of 
Figure |3] where we show 50^ observed sightlines through our 
fiducial disk galaxy. 

The superlinear L^o\-n relation in high critical density trac- 
ers drives a relation between SFR and Lmoi which has index 
less than the Schmidt index (Equations ^Ej . This is simi- 
lar to t he interpretation advocated by Krumh olz & ThompsonI 
(120071) in their models of the SFR-Lmoi relation in GMCs. 

As discussed before, this physical model drives the global 
SFR-Linoi relation in our model galaxies. That said, the con- 
nection between this and how the Ln,oi-n relation is realized 
on a cell by cell basis when considering physical models of 
galaxies, and the effects of radiative transfer and light redis- 
tribution is more subtle. In the remainder of this section, we 
examine the two case examples of n ^ Wcrit of CO (J=3-2) and 
HCN (J=l-0) in detail. 

Non Ground-State Transition — We first consider the case in 
which the gas density is in large part below the critical den- 
sity of the emission line (which is not a ground state transition 
- e.g. CO J=3-2) on a cell by cell basis for our fiducial disk 
galaxy. In the lowest density cells, high critical density trac- 
ers such as CO (J=3-2) are subthermally excited. Unlike the 
case where n ^ rian, the emission from the lower mean den- 
sity cells does not drop monotonically with decreasing mean 
cloud density. Rather, the molecular excitation (and con- 
sequent emission) is supported by radiative excitation from 
neighboring cells with higher mean density. 

The cells with the lowest mean density show a relatively 
constant flux level which is higher than would nominally 
be allowed if collisions alone drove the molecular excitation 
(FigurelDl. To see this, consider the dependence of the cell's 
intensity on the molecular level populations; 

n,iA,ii 
t.= . J, „ , , (8) 

where n„ and «/ are the upper and lower state level popula- 
tions, and Aui, Biu and Z?ui are the Einstein coefficients for 
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Fig. 3. — Velocity integrated molecular line intensities (erg s"' cm"- Hz"' - lim s"') along 50' sightlines for a fiducial disk galaxy versus mean density along 
the line of sight. The crosses represent individual sightlines (with implied spatial resolution of '^ 160 pc) for an individual simulated disk galaxy. The quoted 
mean density is the average along each sightline, and the velocity-integrated intensity is the emission that escapes along each sightline. The disk galaxy has a 
circular velocity of 160 km s"', and a gas fraction of 40% (though the results are general). The solid lines show linearity, and the normalization is arbitrary. The 
CO (J=l-0) emission traces the mean density along the line of sight roughly linearly, while the hig her c ritical density tracers trace mean density superlinearly. 
These results are expected in order to recover the observed SFR-line luminosity relations (Equations l5l7t . 
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Fig. 4. — CO (J=l-0), (J=3-2) and HCN (J=l-0) intensity versus mean cloud density in our fiducial disk galaxy on a cell by cell basis. Each cell typically 
contains a complex of GMCs. The quoted mean density is the mean density of the 50' cells in our sample giid. The blue line shows linearity, and the colors 
in the points show the ratio of LTE level populations to actual level populations for the upper J level of each transition. See text for details on how this drives 
the observed SFR-line luminosity relations. The units of intensity are in erg s" 
luminosity units. 

spontaneous emission, absorption, and stimulated emission, 
respectively. In subthermally populated levels, n; ':^ n^, and 
consequently Ii, oc riu/rii. When both / and u are subthermal, 
«; and n„ increase in lock step superlinearly with increasing 
mean cloud density. This owes to a combination of the effects 
of collisional excitations (which increases linearly with den- 
sity) and the additional contribution of radiative excitations. 
Because «/ and n„ both increase monotonically with gas den- 
sity, there will be little change in the intensity as a function of 
increasing gas density. Consequently, the value of the inten- 
sity is roughly constant at («„Aa;)/(«/B/„) where «„ and n/ are 
inflated (above the effects of collisional excitation alone) by 
line trapping (Figure [Hi. 

At higher mean cloud densities, when the collisions begin 
to contribute significantly to the excitation processes for the 
/ state, n„ continues to increase superlinearly with increas- 
ing density whereas «/ increases only linearly. Thus, «„ in- 
creases with respect to «/ with increasing mean cloud density, 
and the line intensity begins to rise with cloud density (Equa- 
tion[8]l. Emission in these cells is coming both from subther- 
mally excited gas, and the small fraction of thermalized gas 
(in cloud cores) in these relatively low mean density cells. At 
the highest mean cloud densities, the I and u states are both 
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cm - Hz ' which are proportional to the standard observed L' molecular line 

thermally populated, and the emission rises linearly with in- 
creasing mean density as in § 13.2.11 In Figure |5] we plot an 
illustrative example of how the level populations evolve with 
increasing mean cloud density by showing the CO level pop- 
ulation distributions for three cells with different mean densi- 
ties in our fiducial disk galaxy. These rates of level population 
increase with increasing mean gas density translates to the re- 
lationship between Iqo and gas density for the CO (J=3-2) 
transition as shown in Figure |4] (middle panel). 

At first glance. Figure H] suggests that the CO (J=3-2) in- 
tensity traces the gas density (/3) sublinearly, which is con- 
tradictory to what we would expect given the nearly linear 
SFR-CO (J=3-2) index (a; Equations |5]|7]i. However, the im- 
portant quantity to consider is the total integrated intensity 
summed along sightlines through the galaxy (Figure [3]i. The 
gas in the low density subthermally populated regions along 
the line of sight is excited by emission from warmer, higher 
density gas, and thus has a characteristic intensity reflective of 
a higher brightness temperature than the meager densities in 
these gas cells would normally allow for via collisional exci- 
tation alone. Conversely, emission from higher density gas is 
more representative of the dense regions the photons originate 
from. If one could see directly into the thermalized nucleus. 
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Fig. 5. — Normalized level population distributions for three gas cells of 
varying mean density in our sample disk galaxy. In the lowest density cell 
(n ~10~cm"'), all J>1 states are subthermally populated, and rising with 
increasing density. The J=l level is increasing with respect to roughly fixed 
J=0 populations with rising cloud density, resulting in increasing CO (J=l-0) 
flux with rising mean cloud density. Conversely, in the lower density curves, 
while the level populations from higher lying levels are also increasing with 
mean cloud density, the are not increasing with respect to each other (e.g. the 
relatively constant J=3/J=2 level population ratio from the n ~10^cm"^ cell 
through the it '^lO^cm^^ cell). Here, the flux from e.g. the J=3-2 line will 
remain constant with increasing mean cloud density. This occurs until the 
J=2 level approaches thermalization and the J=3 populations begin to increase 
with respect to the J=2 populations. 



the CO (J=3-2) emission would be characteristic of the den- 
sity of gas traced, and the co nsequent /cp-gas density relation 
would be linear (e.g. § 13.2. lb . However, line of sight observa- 
tions of the galaxy contain a contribution both from the dense 
nucleus, and the subthermal cells along the LOS which have 
a rather high L,noi/nceii ratio. The contribution to the emission 
from lower mean density regimes along the line of sight re- 
sults in a total superlinear relation between intensity and line 
of sight mean gas density (f3). 

We show this more quantitatively in Figure |6] where we 
plot the relative light (line flux) to density ratios versus the 
mean cell density along a single sightline peering through the 
nucleus of our sample disk galaxy. We plot the light to den- 
sity ratios for the CO (J=l-0) and CO (J=3-2) transitions, and 
normalize the ratios at the highest density cell. In Figure |6] it 
is evident that the lower density gas along the line of sight to- 
ward the nucleus proportionally emits substantially more CO 
(J=3-2) intensity as a function of cell density than than the 
thermalized gas ^. This translates to the CO (J=3-2) intensity 
tracing the sightline-integrated mean gas density superlinearly 
(Figure |3]l. 

Ground State Transition — We now briefly turn our attention 
toward the SFR-HCN (J=l-0) relation as it represents an in- 
structive combination of the phenomena driving the SFR-CO 
(J=l-0) and SFR-CO (J=3-2) relations. At its root, the charac- 
teristic emission pattern from HCN (J=l-0) falls into a similar 
category as that of CO (J=3-2). That is to say, it is dominated 
by subthermally excited gas, and the emission from these cells 
traces the gas density in a superlinear manner (e.g. Figure |6l 
which was shown for CO J=3-2, but the results are applica- 

' It is important to note that there is no "extra" photon production from 
subthermally excited gas. The emission from this diffuse gas is simply redis- 
tributed light which largely originated in thermalized cores. It is the increased 
relative light to density ratio (Figure|6) in this subthermally excited gas along 
the Une of sight that causes the superlinear /3 index (Equation|7). 
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Fig. 6. — Ratio of intensity to gas density in gas cells along a line of sight 
that goes through the nucleus of the sample galaxy for the CO (J=l-0) and CO 
(J=3-2) transitions. Both the intensity and density curves are normalized to 
their respective maximum values (which both occur in the densest cell), so the 
expected light to density ratio is unity at the maximum density. For the CO 
( J= 1 -0) transition, most of the gas is thermalized along the line of sight, so the 
relative light to density ratio is nearly unity for most clouds. For the CO (J=3- 
2) transition, diffuse gas along the sightline contributes proportionally more 
luminosity than thermalized cells owing to radiative excitations by denser 
gas in the vicinity. When comparing to Figure |2] it is evident that these 
lower density clouds are numerous along a given sightline. The summed 
contribution of the intensity from these low density cells to the total integrated 
intensity causes the relationship between total luminosity and gas density to 
be superlinear. Consequently, the resultant SFR-molecuIar line relation has 
index of order unity. 

ble to HCN J=l-0 as well). Indeed, the result is expected as 
the critical densities from the two lines are only an order of 
magnitude different, and both are substantially higher than the 
mean density of most clouds in our model galaxies. The dif- 
ference between HCN (J=l-0) emission and the CO (J=3-2) 
emission is that there is no flat portion of the line intensity- 
gas density curve as it is a ground state transition (Figure HI 
right panel). The superlinear emission continues monoton- 
ically with cell density across the full range of densities in 
the example galaxy owing to nearly constant J=0 level pop- 
ulations and rising J=l level populations (Equation [8]l. This 
was explicitly seen in Figure |5] for the case of CO (though 
qualitatively the case of HCN is similar). The intensity rises 
superlinearly with increasing gas density because the excita- 
tion processes to the J=l rotational state have a large contri- 
bution from the molecular line radiation field. The source of 
excitations from line trapping increases superlinearly because 
the source of radiation for the excitation - dense, thermalized 
cores - increases superlinearly with n. This superlinear rela- 
tion between Lmoi and n in the emitting gas cells results in a 
nearly linear relationship between the SFR and HCN (J=l-0) 
luminosity in our mode ls, consistent with the we ll established 
observational results of iGao & Solomon! (l2004alfa) . 

We therefore arrive at the following conclusions driving ob- 
served molecular line-SFR relations: 

• For lines with critical densities well below the mean 
density of the clouds in the galaxy, the emission line 
will trace the total molecular content of the galaxy. In 
these cases, we find an SFR-molecular line luminosity 
index equivalent to the Schmidt-law index. This results 
in an SFR-CO (J=l-0) slope of -1.5 when the SFR is 
constrained by SFRoc p'^. Observationally, the SFR- 
CO (J=l-0) index is found to lie between 1.4-1.6. 

• For lines with high critical densities, a superlinear rela- 
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tion exists in Lmoi and n in the gas owing to thermalized 
gas lying only on the high density tail of the density 
distribution. The light from the thermalized cores is re- 
distributed such that the final surface of emission is dif- 
fuse, subthermally excited gas along the line of sight. 
When the line is a ground state transition (e.g. HCN 
J=l-0), the intensity from gas cells rises monotonically 
with mean cloud density, though superlinearly owing 
to heavy contribution to the excitation from line trap- 
ping. When the line is a transition above the ground 
state (e.g. CO J=3-2), emission from the subthermally 
excited cells is roughly constant with increasing mean 
gas density until the level populations involved in the 
transition begin to approach LTE, at which point the in- 
tensity is roughly linear with mean cloud density. In 
either case, for high critical density lines, a superlin- 
ear increase in thermalized gas and high critical den- 
sity photon production with mean gas density results 
in an SFR-line luminosity index (a) lower than the 
Schmidt-law index. In the example of CO (J=3-2) and 
HCN (J=l-0) presented here, this results in SFR-line 
luminosit y indices of ^1, consist ent with the measure- 
ments o f INarayanan et alJ (l2005h and iGao & SolomonI 
(l2004allbh . 

3.3. Implications of Results 

The results presented thus far are a natural solution to the 
observed SFR-CO and SFR-HCN relations. Without any spe- 
cial tuning of parameters, the same simulations are seen to 
additionally reproduce characteristic CO emission line mor- 
phologies, intensities, excitation conditions and effective radii 
of loca l disk galaxies and mergers (Naravananet al. 2006a, 
l2007bl) . as well as quasars at z ^^ 6 (Narayanan et al.ii2007 a). 
More broadly, the same galaxy evolution simulations have 
shown successes in reproducing a large body of characteristic 
observable features of starburst galaxies, ULIRGs and quasars 
from z=0-6 (§ 12.11 and references therein). In this sense, the 
modeled reproduction of the observed SFR-CO and SFR- 
HCN relations are a natural result of a series of simulations 
which have reproduced many observed phenomena. 

Second, we re-emphasize that the relationship between SFR 
and CO/HCN emission is built in to the excitation mecha- 
nisms of both diffuse and dense gas in galaxies, though should 
be taken in an ensemble sense. Because the average physical 
conditions of the disk galaxies and merger snapshots in our 
simulations result in molecular line-gas density relations com- 
patible with the observed SFR-molecular line relations, no 
particular combination of disk galaxies and/or merger snap- 
shots was necessary for the reproduction of observed relations 
in Figure [T] 

A natural question is whether or not any given SFR- 
molecular luminosity relation holds particular significance as 
a physical SFR "law" relating the SFR to a property of the gas 
itself. The results in this section show that all SFR-molecular 
line relations are reflective of the underlying Schmidt law re- 
lating the star formation rate to gas density. That is, our 
models suggest that the observed linear SFR-molecular line 
luminosity relations (for high critical density tracers) do not 
represent a fundamental indicator of the SFR, but are rather 
simply indicative of the underlying Schmidt law. Because 
we assumed that the SFRoc p'^, the resultant SFR-molecular 
line relations from the models matched the observed relations 
rather well (Equations |5]|7]i. To some degree, these models 
suggest that flie observed SFR-CO and SFR-HCN relations 



reflect a physical SFR law similar to the ones assumed for 
this model (with index A^ « 1.5). Moreover, simulations have 
shown that this choice of an SFR law reproduces th e observed 
surface den sity Kennicutt ( 1998a b) SFR laws well dCox et al] 
2006c', Springel i2000l) . In this sense, the predicted molecu- 
lar SFR relations in this section are simply reflective of the 
existing surface density SFR laws, as well as volumetric gas 
density SFR laws. 

Finally, we note that a general consequence of our models 
is that at the highest SFRs, the SFR-Lmoi slope will naturally 
turn toward the underlying Schmidt index. This is because 
at these SFRs, the mean density of the galaxies is typically 
high enough that the mean density is of order the critical den- 
sity of the molecular transition. In these cases, the molecular 
line emission will be essentially counting clouds in a man- 
ner similar to the CO (J=l-0) emission described earlier, and 
an SFR-Lmoi index of ^1.5 will result. This is simi l ar to t he 
interpretation advocated bv lKrumholz & Thompsonl(l2007h . 

4. TESTABLE PREDICTIONS 

In FigurelU we showed for a random drawing of our sample 
of ^100 disk galaxies and merger snapshots that our model 
results displayed consistent SFR- CO and SFR-HC N indices 
with observations of local galaxies (iGao & Solomo n 2004a, 1^; 
iNaravanan et alj|2005l: lYaoet all 12003"). We can further ex- 
tend these model results to make testable predictions for the 
indices of unobserved SFR-CO and SFR-HCN transitions. 

We plot the predicted indices for the range of readily ob- 
servable SFR-CO and SFR-HCN transitions in Figure |I1 To 
derive these results, we randomly drew 35 of our ^^100 galax- 
ies and took the best fitting slope between the SFR and molec- 
ular line luminosity from the entire galaxy. Because we are 
modeling unresolved observations of entire galaxies in our 
simulation sample, these predictions can be directly applied 
to observations. We did this for each transition 100 times, and 
in Figure|2]denote the mean of these results with the solid line, 
and the standard deviation in the dispersion with the hatched 
region. We additionally plot the obs erved data from the sur- 
veys of IGao & SolomonI (l2004"bl) and lNaravanan et alj ( l2005l) 
with associated error bars. The model results are quite con- 
sistent with the three observational data points. Moreover, 
the dispersion in the prediction is consistent with the resultant 
dispersion in slopes of a random sampling of a c omparable 
fraction of galaxies in the lGao & SolomonI (l2004alibh sample. 

The slope of the SFR-CO (J=l-0) relation roughly follows 
the assumed Schmidt law as the CO (J=l-0) serves as an accu- 
rate tracer of the total molecular gas. Higher lying transitions 
become shallower as the relationship between the transition 
and the gas density (globally) becomes superlinear. At the 
highest transitions (e.g. CO > 4), the decrease in the SFR- 
Lco slope begins to flatten. Broadly, this owes to a relative 
plateau in the critical densities as a function of increasing ro- 
tational transition. 

The SFR-L^oi index for tracers with critical density higher 
than HCN (J=l-0) or CO (1=3-2) is directly predicted to be 
sublinear. This owes to the increasingly small fraction of gas 
in galaxies which is thermalized in these transitions. As the 
mean density of galaxies increases, the amount of gas ther- 
malized in these transitions (in the most extreme peaks of the 
density distribution; e.g. Figure |2| increases superlinearly. 
Alternatively said, the fraction of emission dominated by sub- 
thermal emission in galaxies increases as the critical density 
of the tracer increases. Consequently, the Lmoi-n transition 
becomes increasingly superlinear. 
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Fig. 7.— Predicted indices between SFR-CO and SFR-HCN rel ations for a series of mol ecular line transitions. The SF R-CO (J=l-0), SFR-CO (J=3-2) and 
SFR-HCN (J=l-0) indices are all consistent with the observations of Gao & Solomon 1 2O()4aj0) and lNaravanan et aj] 120051) . To simulate observational variance 
with particular galaxy samples, we randomly sampled 35 of our ~100 model galaxies 100 times. The solid line shows the mean derive d slopes, and the shaded 
region the Icr contours in the dispersion of slopes. Additionally plotted are the observed data from the surveys of lGao & SolomonI j2004bl) and lNaravanan et 3] 
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2005) with their quoted error bars. These predicted slopes serve as a direct observable test for these models. 



The predictions seen in Figure |7] may serve as model- 
distinguishing observational tests. That is, the standard in- 
terpretation of the linear relationship seen between Lir and 
HCN (J=l-0) luminosity and CO (J=3-2) luminosity is that 
the molecular line emission traces dense gas more physically 
involved in the star formation process (e.g. n >10^~^cm"^). 
In that picture, observations of HCN and CO transitions with 
even higher critical densities ought to similarly show a lin- 
ear relationship between Lir and molecular line luminos- 
ity as they continue to pr obe dense star forming cores (e.g. 
iGao & SolomonI l2004albl) . In contrast, our models suggest 
that the observed relations between SFR and molecular line 
luminosity are driven by a superlinear increase in Ln,oi with 
increasing n, and will consequently have a sublinear SFR- 
Lmoislope for higher critical density lines than CO (J=3-2) or 
HCN (J=l-0) (Figure|7]i. Observations of alternative CO lines 
(e.g. CO J=2-l) or higher lying CO or HCN transitions in 
galaxies will provide a direct test of these models. 

We note that while the models of Krumh olz & Thompson! 
(12007) do not explicitly predict the SFR-L^oi index for transi- 
tions beyond the ground state, because the physical reasoning 
for their derived SFR-L^oi relation for GMCs is similar to that 
predicted here, in principle similar results can be expected. 
We discuss this work in the context of their models more in 

m 

The dependence of the SFR-Lmoi relations on the relation- 
ship between n and Wcrit lends two other natural tests of the 
models to be borne out. First, galaxies with the highest SFRs 
typically have the largest mean densities. For these galax- 
ies, because n ^ ria-it, an SFR-Lmoi index equal to the under- 
lying Schmidt index (here, 1.5) is to be expected. That is, 
the slope of the SFR-Ln,oi relation should steepen for high 
critical density tracers at the highest SFRs. We remind the 
reader that the SFR-Lmoi relations typically use Lir as a proxy 
for the SFR, which may contain contamination from em- 
bedded AGNs (which we have not included in these mod- 
el s). Thus, salient point of thes e models, as well as those 
of iKrumholz & ThompsonI (|2007|) regarding these high lumi- 
nosity points is that regardless of the potential contribution of 
embedded AGNs these dense, star forming systems will sys- 
tematically lie above the linear relation between Lir and HCN 
(J=1-0)/CO (J=3-2) seen for lower luminosity (Lir< IO'^Lq) 



systems. 

Second, regions within galaxies with higher mean critical 
density will have a different SFR-L^oi index than global mea- 
surements. Lines of sight which probe the dense nucleus of 
an e.g. ULIRG may have densities approaching that of the 
critical density of high critical density tracers. In these cases, 
a superlinear SFR-Lmoi relation may be found. 

5. COMPARISON WITH OBSERVATIONS 

The model results presented here quantitatively reproduce 
the observed SFR-CO (J=l-0), SFR-CO (J=3-2), and SFR- 
HCN (J=l-0) indices for local galaxies. These are found to 
be reflective of the ensemble-averaged excitation conditions 
in these galaxies, and the manner in which the molecular line 
luminosity is related to mean cloud density. 

Surveys of Galactic GMC cores have shown that high 
dipole moment molecules such as CS a nd HCN serve as an 
accurate tracer of dense gas mass (e.g. IShirlev et al.l 120031: 
IWu et al.l l2005h . However, the interpretation from extra- 
galactic sources is mixed. Varying results have come from 
the observational literature regarding the origin of traditional 
tracers of dense gas (such as HCN) fr om extragalactic sur- 
veys. For example, iGreve et aP (l2006h find the HCN emis- 
sion in some systems to arise from dense, thermalized cores, 
whereas Papadopoulos (2007) cite Arp 193 as an example 
in which the bulk of the HCN emission arises from subther- 
mally populated gas. More broadly, constraints on high ly- 
ing transitions i n CO, HCO " *" and HCN from the observa - 
tional samples ofGreve et al.' (2006"), 'Narayanan et al.'('2005K 
[Papadopoulos (2007), Papadopoulos, Isaak & van der Werj 
(120071) and lYao et al.l (l2003h among others evidence a wide 
range of molecular excitation conditions. In this regard, the 
models presented here are consistent with these observational 
results. 

Sensitive observations of high redshift galaxies are be- 
ginning to measure sources in the early Universe in terms 
of t heir place on the SF R-mole cular line relati on as well 
(e.g IGao etal.| 12007: Gre veeFanHo OS: Hainhne et all 120061: 
'Riech ers et al11o06c iWagg et al.n2005^ i2007i: IWeiB et al.1 
2007, for an extensive review, see Solomon & Vanden Bout 
2005 and references therein). Owing to its relatively high 
abundance, most detections at high-z have been in CO (with a 
smattering of HCN, HCO"*" and other molecules), and thus we 
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focus our comparisons with high-z galaxies on CO. A sum- 
mary of a near l y-curr ent sample may be found in Figure 5 of 
lRiechersetai](l2006h . 

In general, the interpretation of SFR-molecular line rela- 
tions at high redshift are muddied by two factors. First, the 
increasing fraction of AGN in high redshift infrared luminous 
sources almost certainly contaminates the Lir from the ob- 
served galaxies and quasars, thus causing a potentially sig- 
nificant overestimate in the SFR. Second, as the objects span 
a large range in redshifts, the rest frame transitions observed 
are typically quite diverse, and thus not always probing the 
same phase of gas. The cumulative Ljr-CO luminosity (over 
numerous trans itions) index fo r sources from 2=0^^6 as re- 
ported by Riec hers et alj (l2006h is ^^1.4, similar to the local 
SFR-CO (J=l-0) relation. Certainly, at lower luminosities 
(LiR<lO'^L0), the molecular line data is dominated by CO 
(J=l-0) observations of local galaxies. The relationship be- 
tween Lir and CO luminosity for these sources then is consis- 
tent with the results of our simulations if the Lir in the lower 
luminosity sources is dominated by dust heating by O and B 
stars. At higher luminosities (Lir> 10'^ L©) the sources are 
characteristically high redshift quasars, submillimeter galax- 
ies and radio galaxies. While the observed CO transition in 
these sources is typically greater than the ground state transi- 
tion, the relationship between Lir and CO luminosity retains 
an index of ^1.4. In contrast, our models (§[3]) would predict 
a slope less than ^1.4-1.6 if a Schmidt law index of 1.4-1.6 
was valid for these high redshift sources. One possible origin 
for this steeper slope is a contribution to the Lir from embed- 
ded AGN. 

Alternatively, the models presen ted here (as well as those 
of iKrumholz & Thompson! 120071) suggest that high critical 
density tracers such as HCN (J=l-0) and high lying CO transi- 
tions may show a natural upturn from the SFR-line luminosity 
relation at the highest luminosities even without the contribu- 
tion of embedded AGN. The galaxies at this end of the lumi- 
nosity range are typically massive starbursts and/or advanced 
mergers with significant amounts of dense gas. As discussed 
previously, and shown in Figure |4] in the high mean gas den- 
sity regime, tracers of dense gas such as CO (J=3-2) and HCN 
(J=l-0) become thermalized. The galaxies in our simulations 
with the highest SFRs (typically recently coalesced mergers) 
contain a large amount of dense gas, and thus most of their e.g. 
CO (J=3-2) emission arises from thermalized cells. The lumi- 
nosity from this gas rises linearly with increasing cloud den- 
sity (/3), and the consequent SFR-molecular line relation (a) 
will be superlinear (nearly equivalent to the Schmidt SFR in- 
dex) for a sample o f these galaxies. This is consisten t with the 
modeling results of IKrumholz & ThompsonI (l2007h who find 
a similar upturn in their SFR-molecular line relations when 
the mean cloud density is much greater than the critical den- 
sit y of the rnolecula r line, as well as the observational results 
of lGaoetal.l(l2007h . 

We reiterate caution, however, that the observed SFR is typ- 
ically derived from Lir measurements, and at these high in- 
frared luminosities, the Lir may h ave a non-negligible con- 
tribut io n from a c e ntral AGN (e.g. Kim, Veilleux & Sanders 
I2002t iTran et all 120011: IVeilleux. Kim & Sanders 2002). 
While in principle our hydrodynamic simulations have the ca- 
pability to investigate the potential contribution from growing 
black holes, a full calculation of the IR SED as well as relat- 
ing it to inferred SFR properties is outside the scope of this 
work, and is deferred to a future study. 



6. RELATIONSHIP TO OTHER MODELS AND INTERPRETATIONS 

We now turn our attention toward other models and inter- 
pretations for the observed SFR-molecular line relations, and 
h ow our solution compares to these works. 

iGao & SolomonI (l2004bl) interpret the tight linear correla- 
tion between Lir and HCN (J=l-0) luminosity as evidence 
for an increasing fraction of dense gas in the most luminous 
sources in their sample, as well as a constant star formation ef- 
ficiency in terms of dense (n >10^cm"-') molecular gas mass. 
Because stars form in the dense cores of clouds, the linear re- 
lationship between Lir (which is powered by star formation) 
and HCN (J=l-0) luminosity is interpreted as a natural one. 

Our models find that the fraction of dense molecular gas 
naturally increases with star formation rate. This is true both 
for the mergers which funnel cold gas into the central kilopar- 
sec (fueling starbursts of ^100 Moyr"'), as well as isolated 
disk galaxies. Indeed this is an expected result as the SFR 
in our simulations is parametrized in terms of the cold gas 
density. This results in a higher fraction of the gas being ther- 
malized in high critical density tracers in the systems with the 
highest SFRs (and, by extension, highest infrared luminosi- 
ties). Indeed, this is what causes the upturn in the SFR-Lmoi 
relation for high critical density tracers at the highe st SFRs. 

UtiH zing observations of Galactic cloud cores, IWu et alj 
(l2005h found that the dense Galactic cloud cores showed 
a linear relationship between Lir and HCN (J=l-0) emis- 
sion. Using these results, they posited that the linear rela- 
tionship between Lir and HCN (J=l-0) emission arises from 
HCN tracing fundamental star forming units which scale self- 
similarly from star forming regions in the Galaxy to ULIRGs. 
In the context that even high critical density tracers such 
as HCN and CS are thermalized in star forming cores (a 
result supported by both observational and theoretical evi- 
dence, e;£^_ShMev_e^ak_2()03]_Walker, Malonev & Seraby^ 
Il994t IWalker. Narayanan & BossI [l994), the model results 
presented here may suggest that a Kennicutt-Schmidt index of 
A^= 1 .5 may not apply to the star formation rate of dense cores. 
The simulations here show that when the mean cloud density 
approaches the critical density of a molecular line, the line lu- 
minosity faithfully traces the cloud density (e.g. Figure|4l left 
panel). Thus, in the context of these models, the relationship 
between the SFR and line luminosity in dense, thermalized 
Galactic cores is expected to take an index of A^=1.5 if this in- 
dex controls the Kennicutt-Schmidt SFR law for these objects. 
It may be, then, that in light of the linear relation between Lir 
and HCN (J=l-0) that a linear SFR law (e.g. SFR^ p) de- 
scribes star formation in dense cores (|Wu et al.l 20051) . 

The recent study bv lKrumholz & ThompsonI (1200 7^ utilized 
radiative transfer modeling coupled with physical models of 
GMCs (consistent with turbulence-regulated star formation) 
to derive a motivation for the observed SFR-molecular line 
relations. These authors found that for individual star forming 
clouds, the line luminosity from high critical density tracers 
such as HCN (J=l-0) increased superlinearly with mean gas 
density owing to a superlinear increase in the fraction of ther- 
malized gas with increasing mean cloud density. This results 
in a linear relationship between SF R and HCN (J=l-0) lu- 
minosity. While the simulations by IKrumholz & ThompsonI 
(12007) were geared toward physical models of GMCs (in 
contrast to the hydrodynamic models of galaxies and galaxy 
mergers presented here), the physical mechanisms driving the 
SFR-Lmoi relation are the same in both cases. In our simu- 
lations, a small fraction of the gas is typically thermalized in 
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high critical density tracers (recall Figure |2]i. As the mean 
density of the distribution of cloud densities increases (e.g. 
as a galaxy achieves a higher SFR, typically through com- 
pression of gas in the nuclear regions), the fraction of gas 
that is thermalized in high critical density tracers moves to- 
ward the mean and increases superlinearly. An observable 
consequence of this in our simulations is the excitation of 
subthermal gas along the LOS, which is readily testable in 
local galaxy samples. In this sense, it is quite appealing 
that two different modeling techniques of systems of differ- 
ent size scales (GMCs versus models for galaxies) find the 
same physical mechanism driving the observed SFR-Lmoi re- 
lation. In both cases, the fundamental SFR relation is indeed 
the Schmidt relation. The model results presented here may 
be directly tested as they predict that tracers of higher critical 
density than HCN (J=l-0) or CO (J=3-2) will be thermahzed 
in a small fraction of the galaxy's gas mass, and show a sublin- 
ear relationship between SFR and line luminosity (Figure|7]i. 
Finally, we note that we include constant molecular abun- 
dances throughout our model galaxies, as a full chemical re- 
action network is both outside the scope of this work and 
not feasible given our spatial resolution limitations. In this 
sense, we are unable to evaluate these models in terms of po- 
tential HCN chemistry in the vicinity of a hard X-ray flux 
as has been argued by some authors, (e. g. lCombesI l2007t 
iGracia-Carpio et al. 2006; Lintott & Viti 200a)- We do note, 
though, that our models quantitatively reproduce the observed 
relations between SFR and CO (J=l-0), CO (J=3-2) and HCN 
(J=l-0) emission in local galaxies while utilizing constant 
fractional molecular abundances. This may imply that poten- 
tial chemistry-related effects have a negligible effect on ob- 
served SFR-molecular line relations. 

7. CONCLUSIONS AND SUMMARY 

We have utilized a combination of 3D non-LTE radiative 
transfer calculations with hydrodynamic simulations of iso- 
lated disk galaxies and galaxy mergers to derive a physical 
model for the observed SFR-molecular line relations. We 
specifically focus on the examples of the SFR-CO (J=l-0), 
CO (J=3-2) and HCN (J=l-0) relations as they are the best 
constrained observationally, and show that our model quan- 
titatively reproduces the observed relations when a Schmidt 
index of '^ 1.5 is assumed. 

While the linear relationship between SFR and high criti- 
cal density tracers such as HCN (J=l-0) and CO (J=3-2) in 
galaxies have been standardly interpreted as a fundamental 
SFR law owing to dense gas being the formation site of mas- 
sive stars, our models suggest that this is not the entire story. 
The linear relations of SFR and HCN (J=l-0) and CO (J=3- 
2) arise as a consequence of a superlinear relation between 
line luminosity and mean gas density in galaxies. This owes 
to small amounts of thermalized gas in high critical density 
tracers. The fundamental relation is instead the underlying 
Schmidt law which sets the way in which observed transitions 
trace the molecular gas. 

Our model makes the prediction that for CO lines with 
Jupper >3 and HCN lines with Jupper >2, the SFR-line luminos- 
ity relationships will be sublinear (Figure |7]i. These predic- 
tions can directly be tested with existing submillimeter-wave 
technology, as well as with ALMA. Our models additionally 
provide specific interpretation regarding the existing observed 
SFR-molecular line relations: 

1. The slope in a given SFR-molecular line luminosity re- 



lation is dependent on both the underlying Schmidt law 
controlling the SFR for the galaxy and the relationship 
between molecular line luminosity and density of emit- 
ting gas (the details of which are outlined in §|3]l. When 
line luminosity traces gas density linearly, the resul- 
tant SFR-line luminosity index is similar to the Schmidt 
law index. In cases where the line luminosity increases 
with gas density superlinearly, the SFR is related to line 
luminosity with an index less than the Schmidt index. 
The relationship between line luminosity and gas den- 
sity depends on how the critical density of the line com- 
pares with the mean density of the bulk of the emitting 
clouds. This is sim ilar to the explanations posited by 
iKrumholz & Thomp son (2007) who utilized models for 
turbulence-regulated GMCs. This directly affects the 
observed relations between SFR and CO (J=l-0), CO 
(J=3-2) and HCN (J=l-0) in the following way: 

(a) Owing to its low critical density, the CO (J=l-0) 
Une is roughly thermalized throughout most re- 
gions of the galaxies in our simulation sample. 
This results in a linear rise in CO (J=l-0) luminos- 
ity with increasing gas density, and a consequent 
SFR-CO (J=l-0) relationship with index similar 
to the Schmidt index. For the case of a Schmidt 
index of ~1.5, our simulations reproduce the ob- 
served relation between SFR and CO (J=l-0) lu- 
minosity. 

(b) The critical densities of CO (J=3=2) and HCN 
(J=l-0) are much higher than the mean density 
of the bulk of the clouds in our simulated galax- 
ies. Because of this, the relationship between Ljnoi 
and n is superlinear, causing an SFR-L^oi rela- 
tion which has index less than the Schmidt index. 
An important observational consequence of this 
is that significant amounts of emission can arise 
from subthermally excited diffuse gas in the vicin- 
ity of denser regions. For the case of a Schmidt in- 
dex of ^^1.5, the observed relations between SFR 
and HCN (J=l-0) and CO (J=3-2) emission are 
recovered. 

2. The emission processes driving the line luminosity- 
density relations (and consequently SFR-line luminos- 
ity relations) are variable such that some galaxies ex- 
hibit mostly thermalized gas for high critical density 
tracers whereas others are largely subthermally excited. 
Generally, the galaxies with higher SFRs have more of 
their gas thermalized which may drive the upturn in the 
LiR-HCN (J=l-0) relations hip for extr e mely h igh lumi- 
nosity sources observed by iGao et al.l (l2007h . This re- 
sult is additionally recovered by models for GMCs by 
iKrumholz & Thompson (2007). 

3. The physical basis for thes e models has a similar 
reasoning as those posited by IKrumholz & ThompsonI 
(12007.) for the SFR-Lmoi relation in turbulence regulated 
GMCs. In this sense, two completely different model- 
ing techniques arrive at similar physical motivations for 
the observed SFR-Lmoi relation in galaxies. 

4. We reemphasize that these models can be directly tested 
via observations of high critical density tracers (see Fig- 
ure |7]i. In particular, these models suggest that obser- 
vations of HCN lines higher than J=l-0 and CO lines 
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higher than CO J=3-2 should bear subhnear relations 
between SFR and Lmoi- 

5. A natural consequence of this model is that at the high- 
est SFRs, the SFR-L,noii"elation will have an index sim- 
ilar to that of the underlying Schmidt index. This is a 
generic feature of our models, and happens regardless 
of the inclusion of AGN. 
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